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ABSTRACT 

We discuss spatially resolved emission line spectroscopy secured for a total sample of 15 grav¬ 
itationally lensed star-forming galaxies at a mean redshift of z ^ 2 based on Keck laser-assisted 
adaptive optics observations undertaken with the recently-improved OSIRIS integral field unit (lEU) 
spectrograph. By exploiting gravitationally lensed sources drawn primarily from the CASSOWARY 
survey, we sample these sub-L* galaxies with source-plane resolutions of a few hundred parsecs en¬ 
suring well-sampled 2-D velocity data and resolved variations in the gas-phase metallicity. Such high 
spatial resolution data offers a critical check on the structural properties of larger samples derived 
with coarser sampling using multiple-IEU instruments. We demonstrate how kinematic complexities 
essential to understanding the maturity of an early star-forming galaxy can often only be revealed 
with better sampled data. Although we include four sources from our earlier work, the present study 
provides a more representative sample unbiased with respect to emission line strength. Contrary to 
earlier suggestions, our data indicates a more diverse range of kinematic and metal gradient behavior 
inconsistent with a simple picture of well-ordered rotation developing concurrently with established 
steep metal gradients in all but merging systems. Comparing our observations with the predictions 
of hydrodynamical simulations suggests that gas and metals have been mixed by outflows or other 
strong feedback processes, flattening the metal gradients in early star-forming galaxies. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: high-redshift — galaxies: 
starburst— galaxies: kinematics 


1. INTRODUCTION 

Statistical surveys of star forming galaxies at the 
peak era of activity, corresponding to a redshift inter¬ 
val z 1.5-3, have been increasingly complemented over 
the past few years by spatially resolved spectroscopic ob¬ 
servations. Early work focused on characterizing global 
trends such as th e evolution of star formatio n rate den¬ 
sity with redshift (|Madau & Dickinson|2014|) , and evolu¬ 
tion of the “main sequence” of star fo rmation rate (SER' 
as a function of galaxy stellar mass (Noeske et al.||2007| 


Peng et al^p010|). Spatially resolved data from in 


egra. 

held unit (ItUj spectrographs has provided com plemen- 


Schreiber et al.||2006t |Stark et al.||2008t |Law et al.||2009t 

Eorster Schreiber et al. 2009| Jones et al. 2010a| Wis- 

nioski et al.||2011t Swinbank et al.||2012| Wisnioski et al. 

20151), the size and 

spatial distribution of giant star form- 

ing regions (Genzel et al. 2008 Jones et al. 

|2010b Liv- 

ermore et al.||2012 

), and radial metal abunc 

lance gradi- 

ents ( 

Cresci et al. 

2010 Jones et al. 2010b 

|Yuan et al. 


2014 Troncoso et al.|2014| fjones et al.|2015|). 2-D kine¬ 

matic oatiTlmd'TEe properties of clumpy star forming 
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regions provide valuable insight into the emergence of 
primitive disks, the role of instabilities, and the gradual 
assembly of central bulges possibly form ed from inwardly 


migrating clumps (Genzel et al. 2011). The properties 


of these clumps provides important evidence of gravita¬ 
tional instabilities in the primitive disks whereas corre¬ 
lations between their ages and radial positions insight 
into the morphological evolut ion of these early systems 


(Eorster Schreiber et al. 12011 ) 


Metallicity gradients also provide an important oppor¬ 
tunity to study the amount of feedback, i.e., the energy 
returned to the surrounding medium by star formation 
and nuclear activity. Energy input can arise from var¬ 
ious mechanisms including supernovae, radiation pres¬ 
sure, and cosmic rays. These effects combine to drive 
large-scale outflows of gas and metals. Outflows are ubiq¬ 
uitously obser ved from systems with high star formation 
rate densities ( Heckman et al.pQQl ), including virtually 


all star for ming galaxies at high redshifts (e.g. [Shapley 


et al.||2003|). Cosmological simulations confirm that out¬ 
flow mass loss rates must be comparable and often larger 
than star formation rates in order to explain measured 
stell ar mass functions and the mass-metallicity r elation 


(e.g. Vogelsberger et al. 
rates of mass and metal 


2014 Crain et al. 2015). The 


OSS are degenerate with rates of 


subsequent accretion and hence are poorly known. How¬ 
ever, this “galactic fountain” cycle may be constrained 
via its imprint on metallicity gradients. To first order, 
high rates of outflow and subsequent accretion will re¬ 
distribute heavy elements, resulting in a flatter gradient. 
Eurthermore, the redistribution of interstellar gas will af¬ 
fect the radial profile of star formation and hence future 
metal production and feedback. Evidence for these ef- 
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fects has been recently observed in the form of flattened 
met al gradients in t he extended disks of nearby galaxies 
(e.g. Bresolin 2011). Several groups have now explored 
how various for ms of feedback affect metal gradients in 
simulations (e.g. Yang fc Krumholz 2012 Gibson et al 


2013[ I Angles-Alcazar et al. ' 20i4[ ). Despite a range of 


methods, there is a clear consensus that stronger feed¬ 
back (that is a higher rate of energy injection from stellar 
winds, super novae, and other sources) results in flatter 
gradients due to mixing of gas and metals over larger 
physical scales. Various prescriptions are able to match 
observations of local galaxies but predict different be¬ 
havior at early times. We therefore seek to constrain the 
degree of gas cycling via feedback, by directly measuring 
gradients at high redshifts. 

Initial IFU survey s targeted modest numbers of galax¬ 
ies beyond 2 ; 2 (Forster Schreiber et al. 2006 Law 


of securing adequate signal to noise per spatial resolu¬ 
tion element. The advent of multi-IFU spectrogr aphs 
(e.g., KMOS on the ESO VLT; |Sharples et al.||2013|) has 
recently led to a signi ficant improvement m survey ca¬ 
pability. For example, Wisnioski et al. (2015) report re¬ 
solved kinematic data for over 100 galaxies within the 
redshift range 2 < ^ < 2.7 using this impressive instru¬ 
ment. However, the multi-IFU approach comes with a 
major limitation in terms of angular resolution. Typical 
076 seeing conditions correspond to a physical scale of 5 
kpc at z 2, whereas the half-light radius of a typical 
L* galaxy at this redshift is only 2 kpc. Clearly only 
the largest and most massive systems can be adequately 
probed with seeing-limited data. In fact, even with adap¬ 
tive optics, the physical resolutio n is only :^1 kpc ([Gerr 


zelet al. 2006 Law et al.l 2009). Jones et al. (2010a 


2013[ ) illustrate the difticulties of correctly interpreting 


velocity fields and abundance gradients from data with 
so few resolution elements across each galaxy. As an ex¬ 
ample, early IFU surveys suggested that compact galax¬ 
ies at z 2 were characterized by dis persion-dominated 


kinematics with little or no rotation (Forster Schreiber 
et al. 2009). Deeper observations with adaptive optics 


showed that most of these sources in fact harbor rotating 
thick disks, with circular velocities commensurate with 
their l ower masses compared to rotation-dominated sys¬ 


tems (Newman et al. 


2013), confirming ea rlier results 


from le nsed galaxies with superior resolution (Jones et al. 


2010a). 

Studying gravitationally lensed systems with adaptive 
optics represents a highly valuable route to addressing 
the physics of galaxy formation, particularly for the less 
massive and more abundant systems at z > 2. Lensing 
magnification enables higher spatial resolution and bet¬ 
ter sa mpling than is othe rwise possible. In a pioneering 
study. Stark et al. (2008) illustrated the potential of the 
lensing approach by securing resolved kinematic data for 
MAGSJ2135-0102, a z=3.075 star-forming galaxy mag¬ 
nified in angular size by a factor of x8 along its major 
axis by a foreground galaxy and galaxy cluster. Using 
the Keck OSIRIS IFU spectrograph with adaptive optics, 
the source-plane velocity field was sampled with a reso¬ 
lution of 120 pc leading to 20 independent points on its 
rotation curve. Subsequently Jones et al. (2010a) stud¬ 
ied a sample of 6 lensed galaxies, demonstrating that a 
high fraction of such < L* systems display well-ordered 


velocity fields, in contradiction to the claims of earlier 
less well-sam pled studies of nn lensed galax ies of similar 
luminos ities (Law et al. 2009). Likewise, Jones et al. 
(|2010b|) derived a radial metal abundance gradient for a 
lensed z = 2.00 galaxy utilizing the spatial variation of 
[N II]/Ha and [O HI]/Ha on 500 pc scales for a system 
with a half-light radius of 2.9 kpc. The necessity of map¬ 
ping weaker lines such as [N H] makes the metal gradient 
aspect pa rticularly challengi ng, but four systems were 
studied in | Jones et al.| (|20I3|) showing a diverse range of 
gradients, poss ibly dependent upo n t he kinematic prop- 
erties (see also Gresci et al. (20I0|)). | Jones et al. (2013) 
introduced a simple model which suggests that metal gra- 
dients and their evolution should provide insight into the 
radial variation in the mass loading factor governing the 
amount of outflowing gas. More recently, a KMOS study 
of a larger sample of 21 galaxies at ^ ^ I suggests the 
metal gradient tends to be less steep in more intensely 
star-forming systems (Stott et al. 2014). In view of the 
diverse results, clearly larger samples are required, par¬ 
ticularly at high redshift where only lensed galaxies pro¬ 
vide the necessary physical resolution. 

This paper presents the results of a significantly en¬ 
larged IFU survey of lensed galaxies. Two practical de¬ 
velopments have motivated u s to ext end the original sam¬ 
ple presented in Jones et al. (2013). Firstly the number 
of lensed targets with known spectroscopic redshifts has 
increased following the CASSOWARY survey of lensed 
star-forming gal axies located prim arily in the Sloan Dig¬ 
ital Sky Survey (Stark et al.||20I3). The second is a sub¬ 
stantial improvement in the performance of the OSIRIS 
IFU spectrograph following its transfer to Keck 1. This 
exploits the more powerful center launch laser which of¬ 
fers a brighter and more compact LGS beacon, ensuring a 
much improved Strehl. The installation of a new grating 
also improved the instrumental throughput by a factor 
of X1.5-2.5 depending on wavelength. 

The aim of the present survey is to exploit recent im¬ 
provements at the Keck Observatory to extend the orig¬ 
inal campaign, obtaining IFU spectroscopy of a larger 
and more representative sample of lensed galaxies. The 
improved throughput has e nabled us to discar d the se¬ 
lection criterion adopted by Jones et al. (2013) whereby 
targets were initially pre-screened with a long-slit spec¬ 
trograph (NIRSPEC) to ensure adequately strong emis¬ 
sion lines, possibly biasing the sample to unusually active 
sources. Our goals are two-fold. Eirst we aim to mea¬ 
sure the radial metallicity gradient on sub-kpc scales for 
a representative sample of 2 : 2 galaxies and to examine 

further the o rigin of the diverse resul ts obtained by var- 


lous groups. Jones et al. (2013) and Yuan et al. (2013) 
argued that much of the discrepanc y ^n ni etal gradient 


measurements, e.g. iQueyrel et al.| (|20I2|) migh t arise 
from poorly-sampled data, whereas btott et al. (2014) 
suggested correlations with the specihc star formation 
may be the cause. Secondly we aim to characterize the 
kinematics and star formation on sub-kpc scales for a 
larger sample. Our well-sampled kinematics and star 
formation morphologies allows us to evaluate the util¬ 
ity and conclusions drawn from compleme ntary larger 
surveys being undertaken, e .g. with KMOS (Stott et al. 
2014 Wisnioski et al. 2015), at coarser ^5 kpc resolu¬ 


tions. 

A plan of the paper follows. In §2 we discuss the sample 
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drawn primarily from the CASSOWARY survey and the 
relevant selection criteria. We present the OSIRIS spec¬ 
troscopic observations and their initial reduction. §3 dis¬ 
cusses the reductions of the resolved spectroscopic data 
into the source plane utilizing the available mass models 
for the foreground lenses. We discuss the kinematic prop¬ 
erties in §4 and the metal gradients in §5. §6 discusses 
the overall results and we summarize our conclusions in 
^ 7 . 


2. DATA 

2 . 1 . Sample and Observations 

During several observing runs in 2013-2014, we ob¬ 
served 11 gravitationally-lensed sources at a mean red- 
shift z 2 usin g the near-infrared integral field spectro¬ 
graph OSIRIS (Larkin et al. 2006) on the Keck 1 tele¬ 
scopes with the laser guide star a daptive optics (LGSAO) 
system (Wizinowich et al. 2006). Prior to 2013, we ob¬ 
served four similar sources using OSIRIS w ith the Keck 
II AO system and presented the results in [Jones et al. 
(2013). Because of the lower system throughput at that 


time, the early four sources were pre-screened with the 
long-slit spectrograph NIRSPEG to have suitably bright 
[N II] and Ha emission lines so as to ensure a reason¬ 
able signal-to-noise with OSIRIS in practical integra¬ 
tion times. After the transfer of OSIRIS to the Keck 
II telescope and the installation of a new grating in 
2012 December, the AO Strehl was improved a nd the 
OSIRIS thro ughput increased by a factor of x 2 (Mieda 


et al. 


2014). Accordingly, to avoid any bias in select 


mg targets, we abandoned the earlier spectroscopic pre¬ 
screening method. The current sample therefore com¬ 
prises 15 sources in total. Ten new sources were selected 
from the GASSOWARY catalog of st a r-form ing lensed 
galaxies as presented in Stark et al. (2013) based on 
their availability during the scheduled observation pe¬ 
riod, their rest-frame UV spectroscopic redshifts in the 
range z = 1.5 — 2.5, and the presence of suitably bright 
proximate tip-tilt guide stars. We sele cted an additional 
cluster-lensed galaxy (Abell 773) fromjBelli et al.|(|2013) 
with similar criteria as the GAS SO WARY sample. The 
eleven n ew sources ex t end t he UV luminosity and SFR 
range in Jones et al.j ([2013 ) with UV absolute magni¬ 
tudes in the range Muv ^ “22 to —18 and SFR in the 
range 1 to 80 M 0 yr“^ (uncorrected for dust). Most of 
the sources are in the sub-L* regime while a few extend 
slightly beyond the L* magnitude limit (Muv < — 21 ). 
The total sample is presented in Table 
Using reddening-corrected star formation rates (SFRs) 
derives from both Ha and ultraviolet continuum mea¬ 
sures (using the methods described below), we can com¬ 
pare the properties of our present sample with those in 
previous studies. Although there is a wide range overall 
(from ^5 to 250 Mq yr“^), the bulk of our sample have 
values in the range 10 to 100 , with a median of ^JOMq 


lensed surveys (Law et al. 

2009 

Forster Schreiber et al. 

2009 

Wisnioski et al. 20 J 

iffDl 

lE). Mear-infrared pho- 


are available for about half the ejects in our sample 


)me( 

(and hence are not listed in Table [R. However, they he 
in the range log :^9.0-9.6, whiS is significantly less 
than for those in the unlensed surveys which generally 


probe systems from 10^^ to 10^^ Mq. Our survey there¬ 
fore samples star-forming g alaxies somewha t abov e the 
main sequence as defined by Behroozi et al. (2013). 

We closely followed the observing technique described 


Jones et al. ( 2010 b). Observations were undertaken in 


one or more near-infrared passbands (J,H, or K) to secure 
spatially-resolved data on the Ha and [N II] emission 
lines. Where practical, we continued to observe H/3^ and 
the [O III] emission lines within a shorter wavelength 
band. The average spectral resolution is R 3600 which 
corresponds to 3 — 6 A across the J, H, K bands. We 
used the 100 mas pixel scale which gives a field of view 
of at least 1.6 x 6.4 arcsec. We took short exposures of 
each tip-tilt star to center the position before moving to 
the lensed galaxy. Each exposure comprised a number 
of 15-minute sub-exposures with a ABAB dither pattern 
ABAB of increment ^2 — 3 arcsec ensuring that the 
target was present in all frames. Exposure times varied 
from 1—4 hours. The seeing during the observations 
varied between 0.4" to 1.5". The median AO-corrected 
seeing was 0.17". 


2 . 2 . Data Reduetion 
We used the late st OSIRIS Data Reduction Pipeline 


)n, f( 


(Larkin et al. 2006|) Hto perform dark subtraction and 


cosmic ray rejection, followed by a direct or scaled sky 
subtraction prior to spectral extraction, wavelength cali¬ 
bration and telluric correction using faint standard stars. 
Adjacent exposures were used as sky reference frames. 
The data cubes from each exposure were finally com¬ 
bined using a cr-clipped mean. 


2.3. Emission Line Fitting 

We fit a Gaussian to the various emission lines to deter¬ 
mine the line flux, velocity, and velocity dispersion. In all 
cases, Ha is the most prominent line and a key indicator 
of the velocity field and star-formation rate distribution. 
A four-parameter Gaussian curve is fit to the Ha line 
for each spatial pixel using a weighted minimization 
procedure. In each data cube, we select a region devoid 
of emission lines and calculate the weight for each wave¬ 
length from a variance over this region, w{X) = U“^(A), 
for use in the Gaussian fitting. We require that each Ha 
detection must be above 5 cr and, where necessary, we 
spatially smooth the data cube with a Gaussian kernel 
of FWHM 3 pixels to increase the signal-to-noise. Given 
a Ha detection, we fit a Gaussian of identical width and 
velocity to the [N II] emission. Where observed, we fit 
the [O III] and H/d emission lines independently from the 
Ha and [N II] emission lines but with a similar procedure. 
The velocities and line widths obtained from the [O III] 
emission lines are consistent with those obtained from 
Ha to within the la uncertainties. The intrinsic velocity 
dispersion is calculated by subtracting the instrumental 
resolution measured from OH sky lines ainst — 60 km 
s~^ in quadrature from the best-fit line width. 


2.4. Flux Calibration, Extinetion and Star Formation 

Rate 

For each galaxy observed in the Ha observation band, 
we use the tip/tilt reference star to calibrate the abso¬ 
lute flux. We fit a PSF to the image of an integrated 


^ Data Reduction Pipeline Version 3.2 
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TABLE 1 
Observation Log 


ID 

z 

Coordinates 

Dates 

Filter 

Lines 


texp 

FWHM 

FWHM 


AHa 

SFR 



RA DEC 

(MM/YY) 




(ks) 

(PSF) 

(Source Plane) 


(mag) 

(Moyr“^) 


This paper 

cswall 

I.4I 

08:00:12 +08:12:07 

3/13,2/14 

Hn2 

Hq, N[II] 

16.2 

0".21 

1.5 X 1.9 kpc 

1.9 


99 + 90^ 

cswaI5 

2.16 

10:09:01 +19:37:23 

2/14 

Kn2 

Hq, N[II] 

9 

0".53“ 

2.2 X 4.4 kpc 

7.6 

0.52 + 0.10 

42 + 5 




2/14 

Hbb 

0[III], H/3 

3.6 






cswaI9 

2.03 

09:00:03 +22:34:08 

3/13,2/14 

Knl 

Ha, N[II] 

14.4 

0".14 

0.4 X 1.2 kpc 

4.3 

0.46 + 0.17 

59 + 11 




3/13 

Hnl 

0[III], H/3 

3.6 






cswa 20 

1.43 

14:41:49 +14:41:22 

02/14 

Hn2 

Ha, N[H] 

3.6 

0".09 

0.1 X 0.3 kpc 

14 

0.27 + 0.34 

6 + 3 




2/14 

Jnl 

Oflll], H/3 

2.7 

0".25 

0.5 X 1.1 kpc 




cswa28 

2.09 

13:43:33 +41:55:13 

2/14 

Knl 

Ha, N[H 


3.6 

0".23 

0.3 X 1.8 kpc 

9.3 


12 + 11 

cswaSI 

1.49 

09:21:25 +18:10:11 

12/14 

Hn3 

Ha, N[ir 


10.8 

0".22 

0.6 X 2.1 kpc 

3.3 


36 + 33 

cswaI28 

2.22 

19:58:35 +59:50:53 

9/13 

Kn2 

Ha, N[ir 


5.4 

0".15 

0.4 X 0.7 kpc 

10 

1.96l0;f, 

250+Ji 




9/13 

Hbb 

0[IH], H/3 

5.4 






cswaI39 

2.54 

08:07:32 +44:10:51 

2/14,12/14 

Kc5 

Ha, N[H] 

10.8 

0".18 

0.7 X 1.2 kpc 

9.7 

o-38i°:f2 

q q+16 




12/14 

Hbb 

0[IH], H/3 

5.4 

0".07 

0.2 X 0.4 kpc 







2/14 

Jn3 

0[H] 


1.8 






cswaI59 

2.30 

22:22:09 +27:45:25 

9/13,12/14 

Kc3 

Ha, N[H] 

7.2 

0".17 

1.2 X 2.7 kpc 

4.6 

0.44 + 0.95 

53 + 47 




12/14 

Hn3 

0[HI], H/3 

3.6 






cswaI65 

2.13 

01:05:20 +01:44:58 

9/13,12/14 

Kn2 

Ha, N[H] 

7.2 

0".16 

0.1 X 0.2 kpc 

42 

i.UC)_o 58 

7+2 

‘ -3 




9/13,12/14 

Hbb 

0[HI], H/3 

5.4 

0".20 

0.1 X 0.2 kpc 




a773 

2.30 

09:17:57 +51:43:31 

3^3, 2/14, 

Kc3 

Ha, N[H] 

12.6 

0".39 

0.1 X 1.2 kpc 

20.3 

2.5 + 0.25 

30 + 7 




12/14 













3/13 

Hn3 

0[HI], H/3 

7.2 

0".41 

0.1 X 1.2 kpc 







Objects published m| 

\Jones et al.\\2013\ 





J0744 

2.21 

07:44:48 +39:27:26 

1/08 

Kn2 

Ha, N[H] 

9 

0 ".ll 

0.3 X 0.8 kpc 

16 

0.53 + 1.2 





2/11 

Hbb 

0[HI], H/3 

3.6 

0".08 

0.3 X 0.7 kpc 




J1038 

2.20 

10:38:42 +48:49:19 

2/11,3/11 

Kn2 

Ha, N[H] 

9 

0".14 

0.4 X 1.6 kpc 

8.4 

0.67+1.30 

38tll 




2/11 

Hbb 

0[HI], H/3 

3.6 

0".14 

0.3 X 1.7 kpc 




JII48 

2.38 

11:48:33 +19:29:59 

2/11 

Kn2 

Ha, N[H] 

9 

0 ".ll 

0.6 X 0.9 kpc 

10.3 

2.94 + 1.05 

210±167 




2/11,3/11 

Hbb 

0[HI], H/3 

3.6 

0".08 

0.6 X 0.9 kpc 




JI206 

2.00 

12:06:02 +51:42:30 

5/10 

Kn2 

Ha, N[H] 

9 

0".18 

0.5 X 3.0 kpc 

13.1 

1.22 + 0.47 

68 V 2 I 




5/10 

Hbb 

0[HI], H/3 

3.6 

0".33 

0.6 X 3.5 kpc 




Note. 

— The 2-D magnification p, 

is defined as 

the ratio between image plane flux and source plane flux. 

uv 

magnitudes and Ha fluxes 


are of the source planes. “ tip/tilt star is a galaxy. ^ For galaxies with no observations in the band, their star formation rates are 
obtained with the weighted mean extinction Ahq; = 0-8 i 1-0 in the sample. 


flux over wavelength of the star and obtain a ID star 
spectrum from the derived PSF. The flux calibration is 
then calculated at Ha wavelength of the galaxy. The 
only exception is CSWA15 whose tip/tilt reference star 
is an extended galaxy and its flux calibration is measured 
from a UKIRT infrared standard star FS26 observed in 
the same night. We crossed check the all Ha flux cali¬ 
brations with the flux calibrations derived from standard 
stars (FS9, FSll, FS19, or FS26) taken at the end of the 
same nights. Flux calibration derived from these stan¬ 
dard stars generally agree to within 25%. 

Ideally we would do the same process to the H^ ob¬ 
servation band observations. However, there are only 
three galaxies with tip/tilt reference stars taken in the 
H/3 bands (CSWA139, CSWA159 and Abell773). Hence, 


we use the flux calibration from the infrared standard 
stars to calibrate the HP fluxes when the information 
from tip/tilt stars are not available. The flux calibra¬ 
tions derived from the three tip/tilt reference stars are 
in agreemen t with those from st andard stars within 15%. 

We follow Jones et~ar] ( 2Q13| ) in the calculation of Ha 
extinction and star formation r ate. In short, we use d the 
Balmer line ratios Ha/HP with Calzetti et al. (2000) red¬ 


dening curve to calculate dust extinction E{B — V). We 
then calculate the extinction of Ha, Ana = 3.33E{B — 
V). The star formation rate is computed from total 
Ha flux corrected for extinction and lensing magniflca - 
tion with the Ha SFR relation from Kennicutt (1998). 
For the three galaxies with no observations in H/3 bands 
(CSWA11,CSWA28, and CSWA31), we use the weighted 
mean and the standard devi ation of Ahq of the other 12 
galaxies in our sample. Ana = 0.8 ± 1.0, to calculate the 
star formation rates. The resulting Ah^ and SFR are 
listed in Table 1. 


3. SOURCE PLANE RECONSTRUCTION 

By observing gravitationally-lensed systems we can se¬ 
cure much higher spatial resolution and sampling for our 
targets than would otherwise be the case. Nonetheless, 
this gain in resolution is only possible by using accu¬ 
rate mass models for the lens system via which our ob¬ 
servations in the image plane can be transferred into 
the (unlensed) source plane. The key to developing ap¬ 
propriate mass models is the correct identification of 
multiply-imaged systems ideally with spectroscopic red- 
shifts. Given the variety of lenses, from SDSS galaxies to 
Abell clusters, surveyed in this study, it has not been pos- 
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sible to adopt a uniform approach to constructing these 
mass models across our sample. 

When Hubble Space Telescope(HST) images or Gem¬ 
ini Science Archive (GSA) images were not available, we 
took photometric BRI images in good seeing with the 
Echellette Spectrograph and Imager (ESI) on the Keck 
II telescopes. In some cases, e.g. GSWA159 , mass mod¬ 
els w ere already available in the literature ( |Dahle et al. 


2013). Depending on the circumstances, three distinct 


methods were used in constructing new mass models as 
detailed below. The methods and imaging data employed 
to develop the mass models are summarized in T able 


We used th e Light-Traces-Mass (LTM) method (Zitrin 


et al. 12015) to develop mass models for GSWAll, 


USWAK7CSWA19, GSWA28, GSWA31, GSWA139, and 
GSWA159. Generally these lenses are dominated by a 
single galaxy with a few associated companions as they 
represent lenses with angular Einstein diameters in ex¬ 
cess of the size of a SDSS fiber (> 3 arc sec). The method 
assumes that the leasing mass distribution of leasing 
galaxies is described by an elliptical power-law, with a 
gaussian core for the brightest central galaxies, scaled 
according to their luminosities. The member galaxies 
are identified spectroscopically and/or via a color cut. 
To represent the dark matter, the mass distribution of 
the galaxies is smoothed by a 2D Gaussian kernel. Mul¬ 
tiple images are then identified from their similar colors 
and/or spectroscopic redshifts, if available, in conjunc¬ 
tion with a preliminary guess for the mass model. In 
many cas es, spectroscopic r e dshifts are available in the 


litera ture (I S tark et al.| |2013] 
et al.|201If and from onr OS 


Dahle et al. 


2013 


Brewer 


RIS observations. In a tew 


cases where redshifts are not available for identified mul¬ 
tiple images, we leave the redshift as a free parameter. 
In constructing the final best-fit mass model, we employ 
a Monte-Garlo Markov Ghain (MGMG) to minimize the 
image position i.e. the total distance between the pre¬ 
dicted and actual position of the multiple images. The 
model consists of six basic parameters: the exponent of 
the power law, the smoothing Gaussian width, the rela¬ 
tive weight between the galaxy component and the dark 
matter component, the overall normalization of the mass 
distribution, and the two external shear parameters. The 
mass-to-light ratio is fixed for all galaxies except that for 
the brightest central galaxies that we left as a free pa¬ 
rameter in each model. 

Some mass models for GASSOWARY sources, 
GSWA20, GSWA128 and GSWA165 were already avail¬ 
able following modelin g developed by M att Auger (MA) 
as part of that survey ( Stark et al.|2Q13 ). In these cases, 
gri imaging from SDSS was used to ht each lensing com¬ 
ponent with a singular isothermal ell i psoid following the 
procedure discussed in Auger et al. (2011). The same 
mass model is fitted to each of the filters simultaneously, 
although the amplitudes of the surface brightness profiles 
can vary. 

To determine how uncertainties in the lensing mod¬ 
els affect our analyses of both the kinematic state of 
each galaxy and its chemical gradient, we empirically 
analyzed a marginal case of the lensing models. Since 
GSWA28 has a relatively high magnification factor and 
contains only one set of multiple images in constrain¬ 
ing the lensing model i.e. the magnification is the least 


TABLE 2 
Mass Models 


ID 

Method 

Image 

Photometric bands 

cswall 

LTM 

Gemini 

gri 

cswal5 

LTM 

Keck ESI 

BRI 

cswal9 

LTM 

HST 

GVIH 

cswa20 

MA 

SDSS 

gri 

cswa28 

LTM 

HST 

GVI 

cswaSl 

LTM 

Gemini 

gri 

cswal28 

MA 

SDSS 

gri 

cswal39 

LTM 

Keck ESI 

VRI 

cswal59 

LTM 

SDSS 

gri 

cswal65 

MA 

SDSS 

gri 

a773 

JR 

HST 

YJH 


well-constrained, we selected GSWA28 as the appropriate 
case. We constructed 20 different mass models following 
the parameters randomly drawn from the MGMG and 
performed the same kinematic properties analysis (sim¬ 
ple disk model fitting) as in Section below. We found 
that in most of the cases, different lens models primar¬ 
ily affect the position of the source relative to the image 
but the morphology and size of the galaxy in the source 
plane is largely the same. Erom the kinematic model 
fitting in each source plane, the derived galaxy inclina¬ 
tions, position angles and centers are in agreement given 
the uncertainties we quote for the best model. The re¬ 
sulting uncertainty in radius at each pixel is < 10%. In 
a few rare cases, the result from kinematic model fittings 
yields a radius that is 30% higher. Eor GSWA28, if we 
have propagated the expected maximum uncertainty of 
10% in radii into the metal gradient calculation, the final 
uncertainty in the derived N 2 ppo 4 metal gradient would 
have increased by only 2% from ±27% to ±29%. Since 
GSWA28 is one of the least certain lens models, we expect 
that the uncertainties due to lensing models are not the 
dominance source of uncertainty and did not propagate 
this uncertainty into the subsequent kinetic modeling. 

Einally, for the rich cluster, Abell 773, the mass model 
is available fr om a detailed study conducted by Johan 
Richard (JR) ([Richard et al. [2011 , 2010 Limousin et al. 
'2012 Livermore et al.||2012 ). In this case, a parametric 


mass model of the central region of the cluster was devel - 


oped using the LENSTOOL package (Julio et al. 


Briefly, we assumed the cluster mass distribu tion 


2007). 


ollows 


a double Ps eudo-Isothermal Elliptical profile (Elfasdottir 
et al.||2007 ) and added one or more central cluster mem- 


bers as smaller scale perturbations to the mass distribu¬ 
tion. As above, optimization is conducted by minimizing 
the predicted and observed positions of the many multi¬ 
ple images in these well-studied clusters. Table 1 sum¬ 
marizes the overall 2-D magnification /i appropriate for 
each source although, in the following analysis, full ac¬ 
count is taken of the spatial dependence in transforming 
our data to the source plane. 


4. KINEMATIC PROPERTIES 
4.1. Methods and Motivation 

The improved spatial sampling made possible by 
studying gravitationally-lensed sources offers the key to 
resolving several uncertainties of interpretation of the 
kinematic and morphological properties of high redshift 
star-forming galaxies. These include discriminating be¬ 
tween single rotating systems and close merging pairs as 
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Fig. 1. — From left to right: source plane Ho; emission map, two-dimensional velocity field, and one-dimensional velocity (red diamonds) 
and dispersion (green squares) for each of our targets. Ellipses represent the source-plane PSFs. Solid lines show the slits used to extract 
the one-dimensional velocity and dispersion with position angles determined from our simple disk model fits (see text in Section ^ for 
details). Black crosses mark the adopted centers of each galaxy. 
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Fig. 1.— continued 
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Fig. 1.— continued 

well as the validity of well-ordered disk models claimed to 
fit a high proportion of objects studied in the literature. 
For this study we also seek to examine possible correla¬ 
tions between the ratio of ordered and random motions 
[yja- typically used as a measure of the extent to which 
a system is or is not dispersion-dominated) with the pres¬ 
ence of a metallicity gradient (see §5) as this might in¬ 
dicate important diagnostics of how disks assemble and 
mature over a key period in cosmic history. For these 
and other applications, clearly determining the reliabil¬ 
ity of the derived values of the rotational velocity and 
the dispersion is crucially important. 

Early LGSAO-based data on lensed sources already 
highlighted the importa nce of securing well-sampled data 
2010a|bD . 


(Jones et al 


Prior to these studies it was 


commonly believed that systemic rotation was confined 
to the galaxies with l arger stellar masses (Genzel et al. 
2006 Law et a'r]|2009 ). However, due to the small 


ber of lensed sources available, lim i ted sta tistics has 
remained a problem. Jones et al. (2010b) presented 
spatially-resolved dynamics for 6 lensed systems claim¬ 
ing 4/6 showed well-ordered velocity fields - a fraction 
consistent with that de termined for larger, mo re lumi¬ 
nous sources. Recently, |Livermore et al. (2015) collated 
the available data on 17 lensed sources including 10 stud¬ 
ied with natural guide star-assisted AO using SINFONI 
on the VLT and 6 with LGSAO and Keck’s OSIRIS, 


mostly from Jones et al.| (2010b 2013), similarly con¬ 
cluding ^60% of these systems with stellar masses in 
the range 4 x 10^ < M/Mq < 6 x 10^^ are consistent 
with rotating disks. This fraction is lower than that de- 
duced from the KMOS^^ survey ( [Wisnioski et al.||2015[ ) 
based on coarser sampled data of a larger (c^30) sample 
of 1.9 < z < 2.7 galaxies where, over their entire sample, 
it is claimed 83% of the sample are rotation-dominated 
with vja >\ and at least 70% can be considered ‘disk- 
like’ systems. 

Our newly-enlarged set of well-sampled resolved spec¬ 
troscopic data provides a valuable test of the above 
claims and we now examine both the validity of a sim¬ 
ple disk model as a fit to our source plane velocity fields 
as well as the extent to which our systems are rotation 
or dispersion-dominated. Our Ha emission maps and 
their associated radial velocities fields are presented in 
Figure [2 While our Ha-based radial velocity fields are 
in many cases indicative of velocity gradients, the asso¬ 
ciated Ha surface brightness distributions have irregu¬ 
lar, asymmetric, and non-disk like shapes. Such irregu¬ 
lar emission line morphologies are common at these red- 
shifts and star-formation rates and contrast with more 


regular distributions in broad-band filters sampling the 
stellar distributions. To verify this, for each galaxy, we 
attempted to extract the total continuum spectrum ad¬ 
jacent to the Ha line. Although the resulting signal to 
noise of the continuum is limited beyond a scale of :^1 
kpc, it provides a valuable indicator of the center of each 
galaxy. 

First we assess the validity of a rotating disk as a 
representation of the observed H a radial velo c ity field 
following the method adopted in Jones et al. (2010b). 
Those authors used an arctangent function as the sini- 
plest model, viz: 


2 R 

V{R) = Vb H —Vc arctan — 

TT Rf 

where R is the radius from the disc center, 

^9 / .xsinO y cosO .0 

= (xcose -ysinef + ( - — - f 

smi 


( 1 ) 


The simple model has seven unknown parameters: the 
inclination i, position angle 0, the position of the disk 
center (ac, Sc), the scale radius an asymptotic veloc¬ 
ity Vc, and an overall systemic velocity Vq. The variables 
X and y represent the distance of each pixel from the disc 
center. The resulting asymptotic velocity is corrected for 
the effect of inclination. 

We use a Markov-chain Monte Carlo (MCMC) method 
to find the best fitting disk model parameters. We fit the 
observed Ha radial velocity map with the model con¬ 
volved with an elliptical point spread function based on 
the reconstructed tip-tilt star image. In the MCMC fit¬ 
ting, we set fiat priors within the parameter ranges as 
follows. The disk center is set to he within the peak of 
the stellar continuum light distribution (which does not 
necessarily coincide with the Ha emission peak). The 
scale radius is confined to lie within the range of the ob¬ 
served Ha map. The asymptotic velocity is < 3 x AV 
where AK is the maximum velocity difference along the 
slit that passes through the disk center for all position 
angles. We allow the full range for inclination and posi¬ 
tion angles i.e. 0° < i < 90° and 0° < 0 < 360°. The 
uncertainty in each derived parameter is calculated from 
its posterior probability distribution. The results from 
these simple disk model fits are summarized in Table 


4.2. Results 

Our data reveals many cases where there are signifi¬ 
cant kinematic deviations from this simple rotating disk 
model. Because our spatial resolution is <500pc, we 
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find the velocity field for the Ha emission line samples 
the local bulk motion of gas and is not smeared over a 
large area as would be the case in coarser-sampled data. 
Typically our resolution elements have velocity uncer¬ 
tainties of only a few kms sec“^ which often results in 
large reduced despite similar r.m.s velocity residuals 
of ('^modei ~ '^data) ^ 30 km s“^ Compared to studies 
sampling ^ kpc scales (e.g. |Swinbank et al.||2012). 

These significant velocity deviations suggest that sev¬ 
eral of our systems are undergoing various stages of a 
merging process. One-dimensional velocity profiles in 
Figure show many galaxies in our sample can be con¬ 
sidered as a perturbed disk where the large-scale velocity 
map exhibits a clear gradient but the velocity dispersion 
peak is offset from our adopted center. In several cases 
(e.g. CSWA19, and CSWA139), two distinct Ha sources 
of different velocities can be discerned, consistent with 
early stages of a merging pair. There is good evidence 
from either the velocity profiles or the velocity deviations 
from our disk model, that CSWA19, 31, 128, 139 and 165 
represent merging systems. 

At this point, it is interesting to conjecture whether 
such complex kinematic patterns would have been dis¬ 
cernible with a lower spatial resolution. As an ex¬ 
periment we re-analyzed two systems, CSWA19 and 
CSWA31, both clearly poor candidates for pure rotation- 
dominated systems, as if they had been observed at lower 
angular resolution. The source plane velocity fields were 
first smoothed with a typical PSF of 0.6", rebinned to 
a coarser resolution of 0.2", con iparable to that achiev - 
able with the KMOS^^ survey (|Wisnioski et al. 2015), 
and re-analyzed with our disk model. The results of this 
experiment are shown in Figure In both cases the 
fit is significantly improved and there is little evidence of 
any departure from a pure disk model. This accentuates 
the importance of optimally-sampled data for a robust 
interpretation of the velocity field. 

The second question posed above relates to the preci¬ 
sion of the derived rotational velocity and, particularly, 
the ratio of systemic to random motions which is such 
an important diagnostic of the maturity of early disk 
systems. To examine this, we first adopt the inclina¬ 
tion and position angle values obtained from the disk 
model fitting. Despite the velocity deviations from the 
disk model discussed above, the posterior probability dis¬ 
tributions of inclination and position angle are generally 
Gaussian in form suggesting reasonably good estimates 
of these parameters. Other methods for determining the 
inclination e.g. analyzing the elliptical shape of the Ha 


or near-infrared light distr ibutions (e.g. Newman et al. 
2013 Wisnioski et al.|2QT^ are less applicable to our less 


luminous systems with irregular forms. As for the posi¬ 
tion angle, our disk model fits are generally consistent 
(within 2a) with values obtained from a ‘peak-to-peak’ 
velocity method where a pseudo-slit placed at the cen¬ 
ter is rotated to find to yield the largest velocity shear. 
With the exception in GSWA165 where the discrepancy 
between the two methods is large and the disk model fit 
is particularly poor, we adopt the position angle from 
the peak-to-peak method {O' in Table |^. 

Our sample can then be split in two subsets. For five 
galaxies (GSWAll, 15, 28, A773, 159) where the disk 
model represents an acceptable fit {Xred ^ ratio 


Vc/a is probably a valid indicator of whether the system 
is dynamically supported. Adopting the standard often- 
used threshold Vc/a > 1 and calculating the intrinsic ve¬ 
locity dispersion a from the weighted velocity dispersions 
of each spaxel along the major axis, all but GSWA 15 
could be considered to be rotationally-supported sys¬ 
tems (see Table [^. However, the Vc/a ratio is clearly 
inappropriate for the remaining galaxies where the disk 
model is a much poorer fit. An alternative statistic some¬ 
times used is the observed velocity ratio ISVj2a ratio 
where AH is the maximum velocity shear obtained from 
the peak-to-peak method. Applying the threshold crite¬ 
rion of AH / 2(7 >0.4 for rotationally-supp orted systems 
adopted by Forster Schreiber et al. (2009) to the entire 
sample, 9 out of 11 of our sources would then be clas¬ 
sified as rotation-dominated systems. Figure shows 
the correspondence between these two measures {vda^ 
AH/2cr) and the goodness of the disk model {xled) 
our sample. Broadly speaking, both above criteria are 
reasonably consistent in selecting rotationally-supported 
systems. However, without the highly-sampled data that 
allows an adequate test of the goodness of fit, a much 
larger number would be incorrectly placed in this cate¬ 
gory. 

We perform the sa me kinematic fitting method to the 4 
galaxies published in Jones et al. (2013) and find consis¬ 
tent categori zation using the xted criterion above. Jones 
et al. (2010a) already performed kinetic model fitting tor 
J0744 and jl 206 and we fi nd very similar results. J1148, 


which Jones et al. (2013) considered as a rotating sys¬ 
tem using a 1-D velocity profile, has a xled ~ 6.5 with 
a very low velocity residual root-mean-square of 7 km/s 
suggesting a well-ordered disk. Finally, J1038 which is 
visibly a merger, has a Xred ^ 380 and is categorized as 
a perturbed disk. 

It is illustrative, therefore, to compare our conclusions 
to the first results of the KMOS^^ Survey which reports 
kinematics res ult of an IFU survey of ^ 200 galaxies at 
0.7 < ^ < 2.7 ( Wisnioski et al.|2015 ). Our galaxy sample 
has an average AH/sim of 148 km s“^ and average ve¬ 
locity dispersion of 60 km s“^, both of which are broadly 
consistent with the range of AH/sim ^ 100 — 200 and 
cr ^ 30 — 80 km s“l in the KMOS^^ data. However, 
according to Figure ^ only 4/11 (^ 36%) of our sam¬ 
ple would be classifiea as rotationally-supported galaxies 
via the combination of both criteria in Figure whereas 
in the KMOS^^74% of z 2 galaxies were classified as 
being rotationally-supported. The higher percentages in 
the KMOS^^ analysis may arise at least in part from 
a combination of the lower spatial resolution where late 
merging systems are mistaken as a regular rotating sys¬ 
tem (Figure as well as the difficulty of fitting a disk 
model in sucn circumstances. 

5. METAL GRADIENTS 

Given the capability of the current integral field instru¬ 
ments, the only practical method to measure metallicity 
in high-redshift galaxies is to use strong-line metallicity 
calibrators. The direct measurement of metallicity (Te 
method) is not feasible for individual galaxies at this red- 
shift since it relies on the measurement of emission lines 
that are too faint to detect. In this work, we use the [N 
H] A6584/Ha ratio (N2) and the ([O HI] A5008/H/3)/([N 
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Fig. 2.— The potential dangers of interpreting kinematic data with poor sampling. Left: Source plane Ha intensity map and velocity field 
in this study for CSWA19 (top row) and CSWA31 (bottom row). Right: Respective simulated Ha intensity map and velocity field indicating 
the deterioration in resolution equivalent to that for non-lensed sources studied with KMOS. In both cases, the complex morphology and 
velocity structure is lost in the poorer sampled data leading to the erroneous conclusion of kinematically well-ordered systems. 


TABLE 3 

Kinematic Properties of the Samples 


Simple-disk model results Peak-to-peak method results 


ID 

^red 

{Hnodel Hlata) 

*(°) 

en 

Vc(km s-i) 

Vela 

Ay (km s“L 

cr(km s ^) 

e'n 

C\Vl2a 

cswall 

17 

40.2 

64ii 

16± 12 



167 ± 53 

94 ±26 

48 

0.89 ±0.19 

cswaI5 

18 

12.9 

25?i 

194 ± 10 

29f^ 


48 ±3 

34 ±5 

205 

0.71 ±0.06 

cswaI9“ 

9.6 

13.6 

4956 

155 ±4 

126if 

i-6fi 

89 ±2 

78± 14 

153 

0.57 ±0.05 

cswaI9^ 

246 

42.3 


156 ±4 

q' 7'7442 

0 11 248 

3 . 62 ;! 

120 ±7 

78± 14 

150 

1.54 ±0.30 

cswa20 

35 

21.8 


162 ±2 

133if 

2.2f:? 

46± 14 

62 ±25 

168 

0.37 ±0.09 

cswa28 

6.6 

25.3 

23il 

129 ± 10 

911^6 

1 fi2-2 

29 ±25 

56± 12 

127 

0.26 ±0.12 

cswa3I 

93 

20.8 

2211 

249 ±3 

199 I 35 

iZZiio 


102 ± 12 

56 ±8 

244 

0.91 ±0.08 

cswaI28 

990 

69.9 

57 ±3 

210 ±2 

5701^8 

8-6Pd 

262 ±8 

66 ±29 

208 

1.98 ±0.44 

cswaI39 

38 

25.6 

24 ±4 

297 ± 15 

265i?o 

4-3?:i 

166 ± 27 

62 ±32 

294 

1.34 ±0.36 

cswaI59 

7.6 

12.7 

1 3^^ 

175 ±9 

141^' 

2.3f:f 

79 ±29 

61 ±22 

164 

0.65 ±0.17 

cswaI65 

I4I 

17.3 

57 ±5 

164 ±5 

o[;:q384 

00^307 

5.91:5 

83 ±28 

44± 19 

218 

0.94 ±0.26 

a773 

17 

53.7 

9937 

ZZ17 

144 ± 12 

176^^^ 

0 q5.3 
0.^2.2 

54 ±22 

45± 10 

158 

0.60 ±0.28 


Note. — “ Results from fitting to the main galaxy only. ^ Results from fitting to both mergers together as one system. 


II] A6584/Ha) ratio ( 03N2) with the calibration s from 
Pettini fc Pa^dMlI - P PQ4, |Steidel gtir] ( [M4| ) - S14, 


and IVIaiolino et al.| ( |2U08[ ) - MUST 


12 + log (0/H) = 8.90 + 0.57 x N2 (N2, PP04) 

12 + log (0/H) = 8.62 + 0.36 x N2 (N2, S14) 

12 + log (0/H) = 8.73 - 0.32 x 03N2 (03N2, PP04) 

12 + log (0/H) = 8.66 - 0.28 x 03N2 (03N2, S14) 

where N2= log([N H]A6584/H(a) and 03N2= 
log(([0 HI] A5008/H/3)/([N H]A6584/Ha)). 

The PP04 relations were primarily calibrated from the 
Te method based on nearby extragalactic HII regions. 


The M08 relations were based on local galaxies that cover 
a wider range of metallicity compared to the PP04 rela¬ 
tions. The S14 relations were based on a subset of the 
same nearby extragalactic HII regions sample in PP04 
relations that are limited to the range of N2 and 03N2 
observed in 2 : ^ 2.3 sample which have shown to lower 
the systematic offset between the two indicators when 
applied to high redshift sample. To counter the possible 
bias of the inferred metallicity due to the higher N/0 
ratios in high redshift galaxies than the ratios in local 
galaxies, we will compare our measured gradients from 
the N2 calibrators to the gradients from the 03N2 indi¬ 
cators which have been shown to b e less dependent on 
the N/0 ratio than the N2 indicator ( Steidel et al.|2Q14 ). 
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TABLE 4 

Metallicity Gradients 


ID 

[N II]/Ha 

Central 12+log(0/H) 
(from N2 PP04) 

N2 PP04 
dex kpc“^ 

N2 S±14 
dex kpc“^ 

N2 M±08 
dex kpc“^ 

O3N2PP04 
dex kpc“^ 

03N2 S±14 
dex kpc“^ 

03N2 M±08 
dex kpc“^ 



Objects classified as rotationally-supported systems (x^ 

< 20) 



cswal 1 

0.14 ±0.01 

8.54 ±0.06 

8.60 ±0.06 

-0.07 ±0.02 
-0.11 ±0.02 

-0.05 ±0.01 
-0.07 ±0.01 

-0.10 ±0.03 
-0.16 ±0.02 




cswaI5 

0.03 ±0.01 

8.16 ±0.01 

-0.04 ±0.01 

-0.02 ±0.01 

-0.05 ±0.02 

-0.02 ±0.01 

-0.02 ±0.01 

-0.01 ±0.01 



8.17 ±0.01 

-0.03 ±0.01 

-0.02 ±0.01 

-0.05 ±0.01 

-0.02 ±0.01 

-0.01 ±0.01 

-0.00 ±0.01 

cswa28 

0.09 ±0.01 

8.10 ±0.01 

0.11 ±0.03 

0.07 ±0.24 

0.17 ±0.59 






8.33 ±0.12 

0.04 ±0.03 

0.02 ±0.01 

0.07 ±0.04 




cswal 59 

0.08 ±0.01 

8.50 ±0.01 

-0.01 ±0.01 

-0.01 ±0.01 

-0.03 ±0.02 

0.00 ±0.01 

0.00 ±0.01 

-0.01 ±0.01 



8.49 ± 0.03 

-0.01 ±0.01 

0.00 ±0.01 

-0.02 ±0.01 

0.00 ±0.01 

0.00 ±0.01 

-0.01 ±0.01 

a773 

0.17 ±0.07 

8.42 ±0.10 

-0.13 ±0.18 

-0.08 ±0.06 

-0.19 ±0.28 

-0.25 ±0.41 

-0.22 ±0.35 

0.01 ±0.30 



8.43 ± 0.05 

-0.05 ±0.02 

-0.03 ±0.01 

-0.07 ±0.03 

-0.03 ±0.03 

-0.02 ±0.02 

-0.08 ±0.02 

Objects with > 20 

cswal 9 

0.09 ±0.01 

8.30 ±0.01 

-0.01 ±0.01 

0.01 ±0.01 

-0.01 ±0.07 

0.12 ±0.04 

0.10 ±0.01 

0.07 ±0.05 



8.26 ±0.01 

0.02 ±0.01 

0.01 ±0.01 

0.02 ±0.01 

0.01 ±0.01 

0.01± 0.01 

0.03 ±0.01 

cswa20 

0.05 ±0.01 

7.89 ± 0.02 

-0.15 ±0.05 

-0.10 ±0.07 

-0.29 ±0.45 

-0.41 ±0.19 

-0.36 ±0.17 

-1.00 ±0.24 



7.97 ±0.04 

0.05 ±0.02 

0.03 ±0.01 

0.09 ± 0.03 

0.03 ±0.02 

0.02 ± 0.02 

0.25 ±0.05 

cswa3I 

0.33 ± 0.06 

8.50 ±0.03 

0.02 ±0.01 

0.01 ±0.01 

0.03 ±0.01 






8.54 ±0.01 

0.00 ±0.01 

0.00 ±0.01 

0.00 ±0.01 




cswal 28 

0.10 ±0.02 

8.57 ±0.03 

-0.10 ±0.02 

-0.06 ±0.01 

-0.12 ±0.02 

-0.09 ± 0.02 

-0.08 ±0.02 

-0.15 ±0.02 



8.51 ±0.01 

-0.04 ±0.01 

-0.02 ±0.01 

-0.05 ±0.01 

-0.04 ±0.01 

-0.03 ±0.01 

-0.06 ±0.01 

cswal 39 

0.17 ±0.08 

8.25 ±0.01 

-0.01 ±0.01 

-0.01 ±0.01 

-0.02 ±0.02 

0.03 ±0.01 

0.03 ±0.01 

0.02 ±0.01 



8.12 ±0.03 

0.04 ±0.01 

0.02 ±0.01 

0.05 ±0.01 

0.02 ±0.01 

0.01 ±0.01 

0.04 ±0.01 

cswal 65 

0.26 ± 0.06 

8.53 ±0.01 

0.03 ±0.01 

0.02 ±0.02 

0.04 ± 0.05 

0.13 ±0.01 

0.11 ±0.01 

0.16 ±0.05 



8.53 ±0.01 

0.01 ±0.02 

0.01 ±0.01 

0.01 ±0.04 

0.01 ±0.02 

0.00 ±0.01 

0.13 ±0.03 


Note. — For each object, the top line contains values derived along each galaxy’s ‘major axis’. The bottom line contains values derived 
from radial binning. 


For each calibrator, we calculate the metallicity gradi¬ 
ent using two methods: from the pixels along the ma¬ 
jor axis and from all pixels binned radially. The galaxy 
center, position angle, inclination, and pseudo slit along 
major axis are the same as those in Section]^ An advan¬ 
tage of the metallicity gradient measured from the pixels 
along the major axis is that it is independent of the in¬ 
clination and thus less model-dependent. However, it is 
based on fewer measurements than the metallicity gradi¬ 
ent measured from all pixels and sometimes not feasible 
when most of the HII regions do not lie along the major 
axis (e.g. CSWA20). The best-fit gradients and central 
metallicities are shown in Table and plots are shown in 
Figure 

In detail, we calculate an average N2 ratio (or 03N2 
ratio) for each radial bin along the major axis (or along 
the annuli) for PP04 and S14 calibrators. We assume 
that metallicity 12 -b log(0/H) is a linear function of 
radius. Hence, the line ratio is an exponential function 
of radius, metallicity gradient, and central metallicity. 
We fit the N2 ratios (or 03N2 ratios) with the weighted 
least-square regression to find the metallicity gradient 
and central metallicity for each galaxy. 

For the M08 calibrator, we calculate the N2 ratio (and 
[O HIJ/H/3 ratio, if possible) for each radial bin and find 
the metallicity for that radial bin with a brute force max¬ 
imum likelihood estimation. We then obtain the metal¬ 
licity gradient and the central metallicity with a linear 
fit. 

We check the consistency between the metallicity gra¬ 
dients derived from pixels along major axis and radial an¬ 
nuli obtained in section[4j We found that the results from 
the two methods are consistent with each other within 
one standard deviation especially when the rotation is 


relatively well-described by a simple disk rotation. The 
deviation is larger when the major axis is short and the 
measurement along the major axis suffers from a poor 
sampling e.g. CSWA20 and Abell773. 

Metallicity gradients derived from the 03N2 calibra¬ 
tors are consistent with the gradients derived from the N2 
calibrators when using radial binning method and com¬ 
paring between the same set of calibrators e.g. PP04 N2 
with PP04 03N2. Although 03N2 calibrators give more 
accurate metallicity measurement than N2 calibrators, 
the consistency in the gradients derived from the two 
calibrators suggests that the metallicity gradients mea¬ 
sured with N2 calibrators should not be heavily affected 
by a possible variation in N/0 ratios in high redshift 
galaxies. Hence, for the three galaxies that we do not 
have information on [O HIJ/H/3, the derived metallicity 
gradients should be relatively reliable. When compared 
between the S14 and PP04 measurements, the metallic- 
ity gradients derived with S14 calibrators are typically 
flatter than those derived with PP04 calibrators. 

In further analysis, we use the metallicity gradients de¬ 
rived from PP04 in radial binning because we can avoid 
the bias from poor sampling in short major axis cases. 
The PP04 calibrator is chosen so that we can easily com¬ 
pared with other observed gradients from previous ob¬ 
servations. 

6. AGN GONTAMINATION AND BIAS FROM LOW SIGNAL 
TO NOISE DATA 

We now examine whether our sample is affected by 
line emission from AGN. Since our metallicity calibra¬ 
tors are only applicable to bounded HII regions ionized 
by hot stars, any emission from nuclear activity would 
give misleading results for central metallicities and gra- 
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Fig. 3.— A comparison of two methods for determining the 
degree of rotational support in our sample plotted against the 
goodness-of-fit in the disk model. Uc/cr from our disk fits and 
Ay/2cr from the peak-to-peak method is considered in the context 
of the reduced Horizontal dashed lines represent the usually 
adopted thresholds, Ay/2(7 > 0.4 and Uc/cr > 1 for rotationally- 
supported systems and the vertical dashed line represents x^ =20, 
below which we consider the disk model to be appropriate. With¬ 
out adequate sampling, a larger fraction of our sample would be 
incorrectly considered to be rotationally-supported. 


dients. Furthermore the presence or absence of AGN is 
interesting in terms of understanding the relationship be¬ 
tween growth of galaxies and their central supermassive 
black holes. Optical spectra of AGN are characterized by 
high ratios of collisional to Balmer lines, and by broad 
line widths arising from outflows and/or broad line region 
kinematics. These signatures have been confirmed in the 
nuclear regions of many gal axies at z 1 — 2 with spa¬ 


tially resolved spectroscopy (Wright et al.|2QlQ| Newman 
et al. 1120141 IForster Schreiber et al.||2Ul4 IGer 


_^ IForster Schreiber et al. l l^Oni ICenzel et al. 

2Q14p . Notably, Genzel et al. ( |2QT4 ) show that the trac- 


tion of galaxies showing signatures of AGN is strongly 
correlated with stellar mass. 

For each of the eight new galaxies in this paper with 
[O IIIJ/H/3 information, we construct a “BPT”:[N II]/Ha 
and [O III]/H;d diagram from the spatially resolved spec¬ 
tra. Most show a distribution of values entirely consis¬ 
tent wit h star-forming regions as seen in other ^ ^ 2 
galaxies (jSteidel et al. 2014 Shapley et al. 2015). We 
thus only show two examples where the interpretation is 
less clear (Figure |^. GSWA15 shows a central region in 
whose strong line ratios significantly lie at an average of 
0.16 ± 0. 06 dex above the m aximum starburst classifica¬ 
tion line (Kewley et al.|2001 ) as shown in the upper panel 
of Figure 3 However, the central [N H]/Ha = 0.03 ±0.01 
and [S H] 1^6716, 6731/Ha = 0.13±0.03 are both low and 
atypical of AGN. These ratios are at the extreme low 


end p robed by current z 2 samples (e.g. Steidel et al. 


2014) and broadly consistent with an extension of the 


star-forming locus seen in these surveys. Alternatively 
the line ratios are consistent with theoretical expecta- 
tions for an AGN w ith gas metallicity of 0.5 solar 
(Groves et al. 2006), but the surrounding (and clearly 
star-forming) regions of CSWA15 indicate significantly 
lower metallicity. If the central line ratios are affected 
by AGN then the gas-phase metallicity gradient would 
be significantly steeper than we infer assuming stellar 
ionization. 

As an additional check for nuclear activity, we fit Ha 
from the central regions of GSWA15 with a double Gaus¬ 
sian profile. The best-fit broad component has a line 
width cr = 153 km s~^ and comprises 29% of the total 
flux (c.f. cr = 59 km s“^ for the dominant narrow com¬ 
ponent). This is similar to stacked spectra of outer disk 
regions in z 2 galaxies (cr 200 km s~^ and ^40% of 
flux in the b road component; |Newman et al.|2Q13 , Genzel 


et al.|2Qi4 ) and is not indicative of AGN-driven outflows 


(cr > 500 km s or broad-line regions (cr > 1000 km 
s“^). The spectrum of GSWA15 is therefore consistent 
with expectations for emission from star formation and 
associated outflows, and does not show strong evidence 
of an AGN. 

The incidence of AGN in our sample is consistent with 
expectations given their stellar masses. We find no clear 
signatures of AGN and only one possible case discussed 
above out of the 12 galaxie s with available BP T diag¬ 
nostics, including those from Jones et al. (2013). An ex¬ 
ample of typical BPT diagrams that are compliant with 
star-forming regions in our sample, CSWA165, is shown 
in the lower panel of Figured The remaining galaxies in 
our sample have low centraljN H]/ Ha and [S H]/Ha ra- 
tios indicating no evidence for AGN. Genzel et al. (2014) 
find similar results using the same niethods: zero secure 
and only two potential AGN in a sample of 17 2 ; 1-2 

galaxies with stellar mass log M^^/Mq = 9.4-10.3, com¬ 
parable to that of our sample. They likewise find an AGN 
fraction < 10% at logM>^/M 0 < 10.3 based on indepen¬ 
dent X-ray, infrared, and radio indicators. Therefore the 
absence of AGN in our sample is not unexpected given 
the stellar mass range probed. 

We then consider the possible bias in measured metal- 
licity gradients due to possibly low signal-to-noise ratio of 
[N H] emission lines. We select GSWA128 to be a fiducial 
sample in this analysis due to its moderately steep gradi¬ 
ent so that we can examine possible degradation of this 
gradient as we increase the noise. We adopt the 2 cr de¬ 
tection limit of the [N H] emission line for CSWA128;the 
limits for other galaxies are similar. We then syntheti¬ 
cally add Gaussian noise to the data cube up to twice 
the detection limit and re-measure the metallicity gradi¬ 
ent with the same procedure. We compare the measured 
gradients when the noise is added to the data cube to 
the gradients measured from the original data cube in 
Figure We find that the gradients measured from the 
noise-added data cubes are consistent with the original 
measurement up to doubling the noise or a quarter of 
the exposure time for CSWA128. Accordingly we can be 
confident that the measured gradients in other galaxies 
are not likely to be affected given our [N H] detection 
limits. 
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CSWA15 



[Nll]/Ha 


CSWA165 



[Nll]/Ha 


Fig. 4.— BPT diagrams for CSWA15 (upper) and CSWA165 
(lower). Gray points are from individual pixels in each galaxy. 
Red points are from pixels whose valu es lie significantly ab ove the 
maximum starburst line (green curve, Kewley et al.| (|2001|)). The 
corresponding region to the red points is sJiown m red m tJie map 
at the cor ner. Blue curve is the best fitted line for the 2 ; ~ 2.3 
samples in|Steidel et al.|(|2014|). Median values in uncertainties are 
shown in tne purple error bars. A few pixels have large error bars 
e.g. the data point at the top of the lower plot. 



£ 0.0 0.5 1.0 1.5 2.0 


Noise Added (x [Nil] detection limit) 


Fig. 5. — Verifying the validity of our derived metal gradients 
for low signal/noise [N II] data. Metallicity gradient residuals are 
shown for CSWA128 as a function of the amount of synthetic noise 
added (see text in Section!^ for technical details). Gradients are 
derived with PP04 method both along the galaxy ’major axis’ (pur¬ 
ple) and using radial binning(green). 


7. DISGUSSION 


Previous observations (Yuan et al. 2011 Jones et al. 


2013 2015) led to a sirnple picture whereby isolated 


galaxies at high redshifts tend to have steep metallic- 
ity gradients which gradually flatten with cosmic time, 
while interacting galaxies have no discernible gradients 



Redshift 

Fig. 6. — The evolution of metallicity gradient with redshift. Red 
squares represent lensed galaxies from the p resent study (LI5). 
Green squares are the four lensed galaxies in | Jones et al.| (|2013l 
J13). Other data point s are measured from flensed slar burst 
dwarf ( [Jones et al. 1^5] J15), a lensed galaxy (I Yuan e t al.||2Qll| 
Yll), non-l ensed galaxies observed with adaptive opti cs dbwmbank 
et al. |20T^ SI2), and Milky Ways planetary nebulae (jMaciei et ai. 
2UU3| 1VIU3). Open and filled symbols represent kinematically well- 
ordered and disturbed systems respectively. Turquo ise lines show 
predic tions from the hydrodynamical simulations of [Gibson et ah] 
([2013[ G13) emphasizing the strong sensitivity to the incorporated 
feedback. 


at all. Conceptually this could be understood if metal 
gradients only became properly established when galax¬ 
ies are kinematically well-ordered, with gradients flatten¬ 
ing as galaxies grow in size. However, the enlarged data 
set presented here shows a large diversity of gradients 
which is incompatible with this picture. The redshift- 
dependent behavior of metallicity gradients in our sample 
is summarized in Figure ^ As in previous work we find 
(negatively) steeper gradients in rotationally supported 
systems compared to the merging and dynamically im¬ 
mature systems. However, the gradients in rotationally 
supported galaxies are much flatte r (> —0.1 dex/kpc) 
than those previously observed in Jones et al. (2013) 
—0.3 dex/kpc), indicating less evolution in the gradi¬ 
ent slope at z < 2. We now seek to explain the diversity 
in observed gradient slopes especially among isolated ro¬ 
tating galaxies. 

One possibility f or the differe n ce see n between this 
work and that of Jones et al. (2013) could be pro¬ 
cedural. In the present sample, we avoided a possi¬ 
ble selection bias towards unusu ally active or metal-rich 
sources. In Jones et al. (2013), integrated [N H] and 
Ha fluxes were hrst measured for a larger sample with 
a slit spectrograph and only those four sources with rel¬ 
atively bright emission lines (particularly [N H]) were 
selected for subsequent study with OSIRIS. Following 
upgrades to the OSIRIS grating and Keck I adaptive op¬ 
tics system, this pre-screening step was avoided in the 
current sample (with the exception of CSWA31). To 
test for this bias, the left panel of Figure correlates 
the metallicity gradient with the integrated [N 
line ratio for both isolated or rotationally-supported and 
kinematically disturbed systems. All galaxies with [N 
H]/H(a < 0.1 have flat gradients (< 0.1 dex kpc“^). 
Moreover, for galaxies with [N H]/H(a >0.1 there is a di¬ 
vergence between isolated or rotationally-supported and 
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merging or dynamically-immature systems with the lat¬ 
ter showing zero gradients irrespective of the metallicity. 
The right panel of Figure [^examines the dependence on 
the integrated sta r formation rate to check whether the 
Jones et al. (2013|) sample is biased to more intense star- 
forming systems; however, no strong trend is revealed. 

A more likely possibility for the diversity of gradients 
at 2; 2-2.5 is variation in feedback strength. Numerical 

simulations suggest that metallicity gradients are highly 
sensitive to fee dback in the form of outflows and “galac¬ 
tic fountains” ([H lkington et al. | 2Q12[ [Gibson et al.j2Q13 


Angles-Alcazar et al.||2U14]), in' the sense that stronger 
feedback (i.e., higher mass loading factors) leads to more 
gas mixing and therefore flatter gradients. This is illus¬ 
trated in Figure which shows evolution of the metal¬ 
licity gradient in a simulated galaxy using two different 
feedback prescriptions. Simple analytical chemical evo¬ 
lution models show that stronge r feedback also res ults in 
lower gas-phase metallicity (e.g. jJones et al.||2013). The 
correlation of gradient slope witn integrated metallicity, 
as shown in the left panel of Figure could therefore 
also be due to a variation in feedback strength. In fact, 
the enhanced feedback scheme in Figure was preferred 
over the normal feedback scheme conside ring the stel¬ 
lar m ass-halo mass relationship matching (jStinson et al. 

While Figure [^provides evidence that gradient slopes 
are affected to different degrees by feedback, the trend 
has a large scatter suggesting that other factors are likely 
important. One obvious possibility is the degree of rota¬ 
tional support compared to random motions which would 
mix the gas and flatten any gradients. This effect is 
clearly seen in the sense that merging and interacting sys¬ 
tems have flatter gradients compared to isolated galaxies. 
Among isolated galaxies, we would naively expect those 
with higher ratios of rotational to random motion (i.e., 
higher u/cr) to have stronger gradients in metallicity and 
perhaps other properties. However, Figure shows very 
little dependence on u/cr. We compare our observational 
results with z = 2 galaxies drawn from the Illustris sim - 
ulation (Vogelsberger et al. 2014 Nelson et al. | |2015r 


a large-volume cosmological hydrodynamical sirnulation 
with baryonic feedback. Metallicity gradients in the sim¬ 
ulation data were constructed based on a linear fit to 
all of the star forming gas within a radius of 10 kpc 
from each galaxys center. The simulated galaxies have 
gas-phase kinematics and negative metallicity gradients 
which are similar to our observations and likewise show 
no clear trend between gradient slope and rotational sup¬ 
port. In essence, galaxies with gas disks supported by ro¬ 
tation have radial metallicity gradients which are similar 
to systems dominated by random motions, in contrast to 
expectations that bulk velocity dispersion should flatten 
any gradients. This suggests that gradients are estab¬ 
lished on timescales of the order of a dynamical time, 
sufficiently short that prominent random motions are not 
effective at mixing the gas. 

To summarize, we find a large diversity of metallicity 
gradients among rotating galaxies at 2; 2. Most rotat¬ 

ing galaxies in our sample have negative gradients while 
interacting systems have flatter gradients, as is found in 


® Illustris data is available through http://www. 
illustris-project.org 


local galaxy samples. Surprisingly we find that the de¬ 
gree of ordered rotation versus random bulk motion has 
no discernible effect on the gradient slope even among 
isolated galaxies, indicating that gradients are formed 
and destroyed on timescales comparable to the galaxy dy¬ 
namical time 10^ — 10^ yr). The steepest observed gra¬ 
dients with slopes ^ —0.3 dex/kpc are less common than 
found in our previous work suggesting that this is a rela¬ 
tively rare or short-lived phase at high redshift. It would 
therefore be interesting to verify the time-dependence of 
steep gradients in the numerical simulations. Finally, we 
find a weak correlation between galaxy-integrated metal¬ 
licity and gradient slope which we interpret as an effect 
of feedback: the same feedback which reduces the gas- 
phase metallicity evidently also flattens the gradients, as 
found in cosmological zoom-in simulations. The diversity 
of metallicity gradients at 2: 2 is therefore likely caused 

in part by di fferences in the rec ent feedback history. As 
we argued in Jones et al. (2013), the evolution of metal¬ 
licity gradients is evidently a sensitive probe of feedback 
in the form of gaseous outflows and galactic fountains. 


8. SUMMARY 

We present spatially resolved kinematics and gas-phase 
metallicity gradient measurements for a total of 15 grav¬ 
itationally lensed star-forming galaxies at 2: ^ 2 based 
on the analysis of strong emission lines of Ha, [N H], 
[O HI], and H/d. 11 new sources were observed with the 
Keck laser-assisted adaptive optics systems and the up¬ 
graded OSIRIS integral field unit spectrograph to which 
we add earl ier data on four sources presented in [Jones 
et al. (2013). With the aid of gravitational leasing, the 
typical source plane spatial resolution for each source is 
< 500 pc which is considerably better than for other sys¬ 
tems studied at similar redshifts. 

We found the following key results: 


• High spatial resolution is crucial in diagnosing the 
kinematic properties and dynamical maturity of 
2: ^ 2 galaxies. We compare our observed veloc¬ 
ity map of interacting pairs of merging galaxies 
with those simulated at the lower resolution typ¬ 
ical in non-lensed surveys and demonstrate how 
easily such sources can be mistaken to repre¬ 
sent rotationally-supported disks in poorer sam¬ 
pled data. Even for galaxies with no clear morpho¬ 
logical sign of interaction, we often find significant 
large deviations in the velocity field compared to 
rotating disk models which would not be apparent 
in seeing-limited data. As a result we observe a 
significantly lower fraction (36%) of rotationally- 
supported systems in our sample than has been 
claimed (:^74%) from larger kinematic surveys un¬ 
dertaken with lower spatial resolution. 


• We find a much higher fraction of 2: 2 galax¬ 

ies have weak or non-existent metallicity gradients 
than in previous studies of smaller samples ob¬ 
served in the same redshift range. It seems unlikely 
that such a change arises as the result of a bias in 
the earlier sample which pre-selected sources with 
stronger [N H] lines and hence increased metallic- 
ity. We observe only a weak correlation between 
the presence of a gradient and the metallicity and 
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total [Nll]/Ha 3 ^^,. Formation Rate (Msun/year) 


Fig. 7. — Correlation between the observed metallicity gradient for the present sample (L15),|Jones et al. (|2013| J13) and|Swinbank et al.| 
(|2012| S12) and the integrated [N II]/Ha ratio (left panel) and the star formation rate (right panel), rilled symbols represent isolated 
galaxies; open symbols represent interacting and/or kinematically-immature galaxies. 


0.4 



AV/2a 


Fig. 8.— The absence of a correlation between the metallicity 
gradient and the degree of rotational support. Data points repre¬ 
sent the present sample and the grey scale represent the output of 
the Illustris simulation for a sa mple of 2 ; ~2 galaxies (|Vogelsberger| 
let al.|20I4| [Nelson et al.|20I5| >. - 

none with the degree of rotational support, the lat¬ 
ter being consistent with the recent predictions of 
the Illustris hydro dynamical simulations. We ar¬ 
gue that variations in gas and metal mixing due 
to feedback most likely play the dominant role in 
modifying metal gradients and thus can explain the 


sizable scatter we see in our enlarged sample. The 
sensitivity of the observed metal gradient to the 
various modes of feedback indicates it will remain 
a promising tool for understanding galaxy assem¬ 
bly. 
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Fig. 9.— From left to right: source plane [N II]/Ho; map, radial metallicity gradient from N2 calibrators, and radial metallicity gradient 
from 03N2 calibrators for subsample of our targets. The gray dots are measurements from each source plane pixels. Green points and lines 
are measurements from radial binning. Purple points and lines are measurements along ‘major axis’. 
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Fig. 9.— continued 
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